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ABSTRACT 

Recent observational evidence for steep dwarf galaxy luminosity functions 
in several rich clusters has led to speculation that their precursors may be the 
source of the majority of gas and metals inferred from intracluster medium 
(ICM) x-ray observations. Their deposition into the ICM is presumed to occur 
through early supernovae-driven winds, the resultant systems reflecting the 
photometric and chemical properties of the low luminosity dwarf spheroidals 
and ellipticals we observe locally. We consider this scenario, utilising a 
self-consistent model for spheroidal photo-chemical evolution and gas ejection 
via galactic superwinds. Insisting that post-wind dwarfs obey the observed 
colour-luminosity-metallicity relations, we conclude that the bulk of the ICM 
gas and metals does not originate within their precursors. 
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1. Introduction 

The existence of a hot, metal-enriched (e.g. for iron, ~ 1/2 solar), gaseous component 
in the intracluster medium (ICM) of clusters of galaxies has been a well-established fact 
for over two decades (see the seminal review of [Sarazin 1986| ). With upwards of a third 
of a given cluster's total gravitational mass locked up in this gas ([White et al. 1993|) , 



understanding the origin of this massive component has been of the utmost importance. 

That the ICM metals are the byproduct of gas that has been processed in galaxies, and 
subsequently ejected, is now widely accepted - indeed, Larson & Dinerstein (1975) predicted 
that observations of metal-enriched gas in the ICM would be a natural consequence of 
Larson (1974) sup ernovae- driven wind model for elliptical galaxies. Galactic winds are 
certainly the current favoured mechanism for ejecting heavy elements (e.g. |Matteucci fc 



Vettolani 198g ; Pavid et al. 199 lj |Ciotti et al. 1991j |Mihara fc Takahara 199j |Matteucci 



Gibson 1995|) , although it should be noted that ram pressure stripping of the enriched 



interstellar media (ISMs) during the initial phase of star formation may also play a role 
(e.g. |Sarazin 1979| , although for persuasive arguments against the "stripping" scenario, the 



reader is directed to |Whfte 199 1| and |Renzini et al. 19"93| ) 



Matteucci & Vettolani (1988) illustrated that by integrating the mass of metals ejected 
from ellipticals of different initial masses over a Schechter (1976) luminosity function 
(assuming some mass-to-luminosity ratio) with faint-end slope a = —1.25, one could 
account for the mass of iron observed in cluster ICMs as a function of cluster richness. One 
drawback to their early work was that it overproduced iron in relation to a-elements, such 
as Mg, Si, and O. This appears to be at odds with x-ray observations of the hot gas which 
show [0/Fe]> +0.2 and [Si/Fe]> +0.15 ( Mushotzky Tggg ). Subsequent work by David et al. 
(1991) and Matteucci & Gibson (1995) alleviates this discrepancy somewhat by allowing 
for a flatter than Salpeter (1955) initial mass function slope, by mass, of x ~ 1. A more 
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exotic mechanism for ejecting the observed mass of metals in their appropriate abundance 
ratios is the bimodal star formation scenario of Arnaud et al. (1992) and Elbaz et al. 
(1995), although a recent photometric analysis of their model Gibson (1996b) highlights its 
shortcomings, at least in its present form. 

While Matteucci & Vettolani (1988) recovered the total ICM iron mass, one of their 
primary conclusions was that the winds could only contribute at most 10 — > 20% to the 
total gas mass - the bulk of the gas was necessarily of primordial origin. 

The past few years have seen substantial advancements made, observationally, in 
determining the nature of the faint end slope of cluster luminosity functions. The gathering 
evidence seems to indicate that while the bright and intermediate luminosity regime is 
consistent with a uniform a of -1.25 to -1.45 (e.g. Ferguson fc Sandage 1991 ), the faint end 



(i.e. Mb ^ — 15) dwarf spheroidal slope (at least in rich clusters) is significantly steeper, 
with values of a ~ —1.8 — > —2.2 being favoured (e.g. priver et al. 19941 ; pe Propris et al. 
1995] Bernstein et al. 1995; pe Propris et al. 1995| ). 



In a most interesting recent paper, Trentham (1994) speculates that the precursors to 
these dwarfs which currently populate the steep faint-end luminosity functions may in fact 
be the originating source for ~ 100% of the x-ray ICM gas. Trentham (1994) argues that 
if the faint-end dwarf slope is indeed a ~ —1.9, and the precursors eject ~ 8 — > 33% of 
their initial total (i.e. "initial total" equals initial gas plus dark matter) mass in the form 
of gas (i.e. 7 ~ 0.08 — > 0.33), after an initial burst of star formation, then these "dwarf 
precursors" can account for most of the present-day ICM gas, at odds with the earlier 
studies which favoured a primordial origin. 

Following Trentham's (1994) suggestion, Nath & Chiba (1995) derived analytical 
expressions for the amount of mass which can be lost by dwarf galaxies as a function of 
galactic mass. They showed that the Trentham (1994) idea was difficult to support if their 
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more realistic 7 values were considered, but they also did not adopt self-consistent models 
of galaxy evolution. 

In light of the aforementioned observations, and subsequent work of Trentham (1994) 
and Nath & Chiba (1995), we plan to re-examine the potential role played by dwarf galaxies 
in the enrichment of gas and metals in the ICM of galaxy clusters. To this end, we have 
constructed a self-consistent, coupled photo-chemical evolution package suitable for the 
study of spheroidal star systems in the context of a supernovae-driven wind framework. 
This is the first time, to our knowledge, that ICM enrichment models have been generated 
in conjunction with the photometric evolution of the underlying galactic population. Details 
regarding the mechanics of the population synthesis and chemical evolution implementation 
are given in Gibson (1995a,b). Earlier versions of the code were demonstrated in Gibson 
(1994a,b). 

We begin Section ^| with a review of the observational constraints for the problem at 
hand. Following this, we shall introduce our favoured "template spheroidal models" , which 
in turn form the basis for the subsequent analysis of the ICM abundances and gas mass. 
The discussion of our results, and summary, can be found in Sections 3 and 4, respectively. 



2. Analysis 



2.1. Observational constraints 



Previous models which use supernovae-driven winds from ellipticals to enrich the 
ICM have had at least one major drawback - no self-consistency check on the implied 
photometric properties of resultant galaxies (e.g. [Matteucci fc Vettolani 1988| ; |Pastor et al. 



[T989t |David et aTTggg ; |Ciotti et al.T991| ; |Arnaud et al. 1992| ; |Okazaki et al. 1993| ; |Mihara| 
fc Takahara 1994] Elbaz et al. 1995| ; [Matteucci fc Gibson 1995|) . It is all very well to adjust 
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various input ingredients to the models in order to maximise the mass of ejecta and/or 
favour specific abundance ratios, but it is imperative to test that this has not been done at 
the expense of replicating the observed colour-luminosity-metallicity relationships. 

Figure [I] shows the metallicity versus V-band luminosity relation for dwarf through 
giant ellipticals. Our model curves will be discussed in Section |2]^. The absolute magnitudes 
for the Virgo cluster ellipticals were derived assuming H = 85 km/s/Mpc, and indeed 
this value of the Hubble constant is assumed throughout the remainder of the paper. The 
metallicities from the Terlevich et al. (1981) and Sil'chenko (1994) samples were derived 
from Mg2 line index measurements, whereas the lower mass dwarf spheroidals (taken 
from Smith 19~85] compilation) are typically estimated from giant branch locations in the 



spheroid's colour-magnitude diagram (fSmith 1985 , and references therein). As an aside, the 



dominant elemental components of "Z" for these dwarfs are the a-elements, whereas both 
iron and a-elements contribute to the more massive systems flMatteucci 1994|) . Ideally one 
would, for example, like to compare synthetic Mg2 indices for our model giant ellipticals 
against those shown in the figure, but as this is outside the scope of our current analysis, 
we settle for a compromise comparison with the global metallicity Z. In order to estimate 
the value of Z for those galaxies whose abundances were determined with the Mg2 line, (i.e. 
all those points which lie above [Z]=-0.6), we adopt the [Mg/Fe]=+0.25, Mg 2 -Z calibration 
of Barbuy (1994). In light of the measurements of Mg overabundances relative to Fe, of this 
order, in giant ellipticals ([Worthey et al. 1992| ), we felt that this was justified. Our results 
do not change substantially if we were to adopt the older [Mg/Fe]=+0.00 calibration. 



EDITOR: PLACE FIGURE |T| HERE. 



Figure |2| shows the second of our primary observational constraints - the colour- 
luminosity relation. We will be restricting the discussion which follows to the 
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optical-infrared V-K colour, although it should be understood that others, including B-V, 
have been considered. The reason we concentrate on V-K is that it provides a more valuable 
constraint than say B-V. Over the range of luminosities considered in our study, B-V does 
not vary by more than ~ 0.2 magnitudes, whereas the V-K versus My relation is flatter 
and spans > 1 magnitude. The near vertical distribution of ellipticals in the B-V versus 
My plane (e.g. Figure 7b of |Arimoto fc Yoshii 1987]) makes this colour, by itself, a poor 



constraint. 

The sample of giant ellipticals in Figure |2| is taken from Bower et al. (1992), 
and, following their prescription, we shift the Coma data to the Virgo scale using 
A(V — My) = 3.58 and Hq = 85 km/s/Mpc. The local dwarf colours are from Thuan 
(1985). Using the local dwarfs as a constraint can be dangerous, as it is apparent that some 
have suffered complicated star formation histories (e.g. Phoenix, Leo I, Fornax, Carina ... 
Ferguson fc Binggeli 1994] ), whereas others (e.g. Sextans, Sculptor, Draco, Ursa Minor) 



show little or no signs of star formation besides the initial burst. We feel reasonably safe 
in ensuring our model colours trace the lower envelope of the distribution, and simply note 
that subsequent strong star formation epochs will scatter the colours redward in this figure 
(due to the fact that the later bursts will be taking place in pre-enriched material from the 
dying stars in previous bursts, so despite being younger, their colours will almost certainly 
be redder ... see, for example, Table 5A of |Worthey 199 



EDITOR: PLACE FIGURE | HERE. 



We mention in passing one final constraint pertaining to the properties of the elliptical 
galaxies themselves. As Worhey et al. (1992) have shown, and alluded to earlier, giant 
ellipticals seem to possess a magnesium to iron overabundance, as compared to the solar 
ratio, with values between [Mg/Fe]~ +0.2 — > +0.3, albeit with a large scatter. This is 
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almost certainly due to the chemical enrichment history of said systems being dominated 
by Type II SNe, as opposed to Type la. 

Another clue to the importance of Type II SNe, not only for the role they play in 
driving the stellar [Mg/Fe] to super-solar values, but also for the one they play in the 
ICM itself, comes from x-ray observations of a-element abundances in the hot gas. Early 
observational work only allowed detection of the strongest Fe lines ( [Mitchell et al. 1976| ), but 
it has become clear over the past few years, and especially with the results from the ASCA 
satellite, that there is also an a-element overabundance with respect to iron, compared 
with solar ratios, in the ICM gas, with [0/Fe]« +0.2 -> +0.6 and [Si/Fe]« +0.1 -> +0.5 
QMushotzky 1994"|) Again, these ratios are indicative of Type II SNe dominated origin, as 



opposed to SNe Type la. 

The combination of the a-element to iron ratios in both the ICM and the cluster 
ellipticals both strongly implicate early Type Il-driven winds. 

As to the absolute masses of both iron and gas implied by the x-ray observations, we 
refer attention to Figures ^ and^. Arnaud et al. (1992) have demonstrated that a correlation 
exists between the total luminosity originating in a given cluster's elliptical+lenticular 
population and the measured mass of iron residing in the hot ICM gas. The shaded region 
of Figure ^ encompasses the observed scatter in the iron-cluster luminosity relation, as 
reported in Arnaud (1994). Again, a reasonable amount of scatter exists, but the trend 
does appear to be real. For comparison, we note the location of Virgo and Abell 2199, 
poor and rich clusters, respectively. Figure |4| shows the parallel correlation, but this time 
for ICM gas mass, as opposed to iron mass. The gas masses are claimed to be accurate to 
within a factor of two by Arnaud et al. (1992), and that is reflected in the width of the 
scatter at a given luminosity. 
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EDITOR: PLACE FIGURE | HERE. 

EDITOR: PLACE FIGURE | HERE. 

In summary, the primary observational constraints that we must honour in a study of 
this nature are: the colour-luminosity-metallicity relations highlighted in Figures [l] and |2|, 
the [Mg/Fe]~ +0.2 — > +0.3 overabundance in the stellar populations of giant ellipticals, the 
correlation between ICM iron mass and cluster luminosity as seen in Figure and finally, 
the [0/Fe]~ +0.2 — > +0.6 overabundance in the ICM. As this is the first study to couple 
ICM abundances with the photo-chemical evolution of the underlying stellar population, we 
will demonstrate that previous work has suffered due to their restricting of constraints to 
the chemical properties alone. 

2.2. Photo-chemical evolution models 

We utilise the Metallicity Evolution with Galactic Winds chemical evolution package, 
MEGaW flGibson 1996a| ), which is similar in spirit to that of Matteucci & Tornambe 
(1987), but adopts the more aesthetically pleasing (at least in the opinion of the first 
author!) "mass in/mass out" formalism, similar to that of Timmes et al. (1995), as opposed 
to the matrix form of Talbot & Arnett (1973). Diffuse dark matter halos, and their influence 
upon the system's global and gaseous gravitational binding energy, are included, following 
Bertin et al. (1992). As in Matteucci (1992), we adopt initial dark-to-luminous masses, 
and radial extents, of 10. This means that the dark matter halos, although heavy, are very 
diffused and their effect on the potential well is not important, as shown by Matteucci 
(1992). Therefore, we are in a situation of almost minimal binding energy for the galaxies. 
As a consequence, the amounts of matter restored by galaxies into the ICM are close to the 
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maximum ones. Moreover, the assumed constancy of the ratio between dark and luminous 
matter from galaxy to galaxy is leading to a situation where the dwarfs can contribute to 
a maximum amount of matter. In fact, there are suggestions flKormendy 1990| ) that the 
percentage of dark matter in dwarf is likely to be higher than in giant galaxies. We do not 
make any assumption concerning the nature of dark matter which could be either baryonic, 
non-baryonic or a mixture of both. In any case, the nature of dark matter is not relevant 
to the results of this paper. 

A Schmidt (1959) star formation law (i.e. one in which ip varies with some power of 
the gas mass) of the form 

m = umi (i) 

with k = 1, is assumed for the pre-wind (i.e. t < tew) phase, whereas ip{t) = for 
t > tew- The star formation time-scale v is used as a free parameter in order to ensure the 
colour-metallicity-luminosity relations are recovered, as will be shown in Section 3. This is 
similar to the procedure followed by Arimoto & Yoshii (1987) and Yoshii & Arimoto (1987). 
Generally though, v is found to increase with decreasing mass in a manner reminiscent 
of models whose initial time-scale for star formation is set by the mean collision time of 
star-forming fragments in the proto-galaxy ( Arimoto fc Yoshii 1987 ). 



MEGaW has the flexibility to use any number of input ingredient sources. For the 
purposes of this work, we have chosen a universal stellar initial mass function (IMF) of 
slope x = 0.95, consistent with the value implicated by our earlier work ( Matteucci fc| 
Gibson 1995|) , with corresponding lower and upper limits of 0.2 M ( Paresce et al. 1995| ) 



and 65 M , respectively. The main sequence lifetimes come from Schaller et al. (1992), and 
the remnant masses are based upon the analytical expressions of Prantzos et al. (1993). 

For the nucleosythesis yields, we use the most recent metallicity-dependent tables of 
Woosley & Weaver (1995) for masses m > 10 M (i.e. Type II SNe). The classic models of 
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Renzini & Voli (1981) are adopted for single low and intermediate mass stars (i.e. m < 8 
M Q ). Following Iwamoto et al. (1994), we assume that stars in the initial mass range 
8 < m < 10 M undergo core collapse, as opposed to thermonuclear explosion, thereby 
trapping their newly synthesised metals in the resultant remnant (i.e. they only enrich the 
ISM via pre-collapse stellar winds). The binary Type la SNe model of Whelan & Iben 
(1973) has been included, following Greggio & Renzini (1983). For their yields, we use the 
updated Model W7 of Thielemann et al. (1993). We have ensured that the present-day rate 
of Type la SNe is consistent with the observed rate in giant ellipticals (Rj a = 0.07 — > 0.23 
SNu, for Hq = 85 km/s/Mpc [Turatto et al. 1994]) , by choosing, a posteriori, that ~ 3% 
of the mass in the range 3 — ^16 M of the IMF gets locked into Type la-progenitor 
binary systems flGreggio fc Renzini 1983|) . An average value for i?i a of ~ 0.12 SNu for the 



present-day rate in our giant elliptical models was found. 

For gas to be expelled from a galaxy we require the thermal energy of the gas heated by 
supernovae and stellar winds to exceed its gravitational binding energy ( [Larson 1974|) . The 
stellar wind energy, while not important for giant ellipticals, can contribute non-negligibly 
for low mass (i.e. M < 10 9 M ) systems ( Gibson 1994a|) . We use the energy formalism 
outlined in that paper for its inclusion. 

Supernovae remnant interior thermal energy (i.e. the energy assumed to be available 
for driving the wind) evolution follows that of Model B 3 of Gibson (1994b), which draws 
heavily upon the work Cioffi et al. (1988). Unlike the older Cox (1972) and Chevalier 
(1974) equations, which have been used exclusively to date, Cioffi et al. (1988) include a 
sophisticated treatment of radiative cooling in the SNR interior, as well as the effects of 
metallicity upon their evolution. Also, virtually all previous galactic wind models have 
assumed that SNRs evolve in isolation (see Model A of |Gibson 1994b| - in particular, that 



after the initial adiabatic expansion phase a radiative cooling phase is entered and the 
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energy e cools as e(t) oc t a62 , ad infinitum. This is clearly not the shells either 

come into pressure equilibrium with the local ISM, thereby halting the expansion cooling 
term, or, more importantly, come into contact and overlap with neighbouring expanding 
shells, further cooling thereafter being negligible. This was noted in Larson's (1974) original 
paper, but the "evolution in isolation" formalism still stood (and still does in most wind 
models) until the recent work of Babul & Rees (1992) and Gibson (1994b). As we shall 
see, the wind epochs favoured by Model B 3 are significantly earlier than those predicted by 
Model Ao. Again, further details can be found in Gibson (1994b, 1995a). 

The post-tow energetics situation is somewhat more problematic. Scenarios in which 
continuous Type la-driven SNe winds ensue until the present, temporarily driven winds, or 
even no post-tow wind whatsoever, are all feasible. Of great significance in determining 
this outcome is the amount of residual thermal energy which is left in the system after the 
bulk of the gas has been ejected in the global wind at tow- This has been demonstrated 
quite graphically by Arimoto (1989) and Ferrini & Poggianti (1993). We do not wish to 
belabour this point, and as such, we present three different scenarios for the post-tow 
thermal evolution. A minimal model in which the only gas ejected is that due to the global 
expulsion at tow; a maximal model in which all post-tow ejecta from lower mass stars is 
continuously swept out of the system due to the continued heating from Type la SNe; and, 
a standard model, which incorporates Arimoto's (1989) assumption that the fraction of 
residual thermal energy remaining after tow compared with that before is ^ 0.01. These 
latter models usually lead to temporary winds of duration < 0.5 Gyr for massive systems, 
while for low mass systems they coincide with the maximal models. 

In order to calculate the coupled photometric properties of our spheroids, we have 
developed a simple population synthesis package ( pibson 1996b| ) which, for the work 
described here, used the metallicity-sensitive isochrones of Worthey (1994). This compilation 
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spans ~2.5 dex in [Z], and covers the primary evolutionary stages from the zero age main 
sequence to the onset of the post-asymptotic giant branch (or carbon ignition, depending 
upon the initial mass in question). Luminosity weighted (V-band) metallicities [< Z >]v, 
[Mg/Fe], colours and luminosities, and final mass-to-luminosity ratios were all generated as 
described in Gibson (1996b). 



2.3. Intracluster medium implications 

The work of Arnaud et al. (1992) has shown that there is a direct correlation between 
a cluster's elliptical population luminosity, and the measured abundance of iron in the ICM. 
With knowledge of a given spheroid's ejected mass of gas (or element i) and its V-band 
luminosity, we can easily integrate over the cluster's Schechter (1976) LF, as opposed to 
working with the mass function and assuming a priori some typical mass-to-luminosity 

■p /on 

ratio. By normalising to a cluster's E+SO V-band luminosity L v , we can then write the 
total mass of element i ejected into the ICM by galaxies of luminosity greater then Ly in as 

* j \T * 



Mf =LT J >*.(^)V^vd(^)/ 

r min / r * ' 
n V ' L V 

°° T T 



r mm / r * 

V ' V 

where rrii is the mass of element i ejected by an elliptical of luminosity Ly 1 . a is the 
usual faint-end slope, with canonical values ranging from -1.00 to -1.45 (e.g. [Ferguson fc 



3 If the function m>i can be written in the form of a power law in luminosity, then equation 



can be written as a simple incomplete gamma function (e.g. Matteucci fc Vettolani 1988 



Elbaz et al. 1995] ). This is sufficient provided one restricts the fit to the high mass end of 



the LF (e.g. M g (0) > 10 9 M ). This simple power law form does not extend down to the 
dwarf (i.e. M g (0) ^ 10 8 M ) regime. If we were to fit a power law to the high mass end 
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Sandage 1991, and references therein). L v , for the bright-end of the LF, is chosen to be 
L*y = 1.54 x 10 10 L (i.e. M* = -20.6), again, typical for H = 85 km/s/Mpc. 

The recent years have led to a revolution of sorts in our picture of the faint-end 
of the cluster and field LF, and in particular, the realisation that the faint-end is not 
well- described by a single slope a. The ten low and medium redshift clusters in the survey 
of Phillips et al. (1996) show a surprising uniformity in that the faint-end slope brighter 
than My ~ —17.3 is consistent with a ~ — 1, whereas fainter than this, a clear up-turn, 
with a power law slope a ~ —1.5 — > —1.8, is seen. The deep, very faint, LFs, from the cores 
of four rich clusters, presented by De Propris et al. (1995), are even more extreme with a 
faint-end slope a ~ —2.2. The models of Babul & Ferguson (1996) also predict a similarly 
steep slope below My fh —17.3. In Section 3, we consider a number of LFs including 
both the canonical single-slope faint-end, and, what appears to be, the more appropriate 
two-component form, consistent with the aforementioned observations. 

One last point which should be made regarding equation is the adopted lower limit 
on the intergals, My m . In general, or at least for all models in which the faint-end LF slope 
is a ^ —1.45, the results of Section 3 do not depend sensitively upon My m , and thus we 
typically assume the luminosity associated with the lowest mass dwarf in our study (i.e. 
M g (0) = 10 4 M ). As we shall see, though, for some scenarios in which the faint-end slope 

of our standard models in Section 3, and then blindly extrapolate to the low mass end, we 
would overestimate the predicted masses of ejected gas and metals from these dwarfs by 
factors of two to three. Indeed, the predicted ejected gas mass would actually exceed the 
mass of gas initially present! For "normal" a > —1.45 LFs, this unphysical situation may 
not be numerically significant, but as we are interested in exploring the relevance of very 
steep (a « — 1.9 — ► —2.2) faint-end slopes, we feel it prudent to integrate over the Schechter 
function properly. 
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of the LF is very steep, we have to resort to more subtle means in order to set My 111 . The 
reason for doing so is that as Melnick et al. (1977) and Thuan & Kormendy (1977) have 
demonstrated, the maximum fraction of a cluster's luminosity that is tied-up in the LF, 
below My « —17, is ~l/4, based upon cluster diffuse-light constraints. Occasionally, a 
lower integration limit in excess of the "default" must be enforced in order to ensure that 
we do not exceed this ~l/4 dwarf luminosity fraction. At some level, this should not be 
surprising, as one need only refer back to equation (22) of Schechter's (1976) seminal paper 
to see that for a < —2, a regime we are exploring in this paper, the integrated cluster 
luminosity diverges. 

3. Discussion 

We now present the results for our so-called "template models" , the input ingredients 
for which are described in Section [27|. Table |l] lists the output for the three post-tcw 
wind scenarios: the extrema (continuous Type la SNe-driven winds to the present and 
suppression of all la-driven winds) and the standard model (temporary winds which die at 
after ~ 0.2 — > 0.5 Gyr). For each model we list - column (1): the initial luminous mass in 
gas (in M Q ) - column (2): the star formation astration parameter of equation [1] (in Gyr -1 

- column (3): the time of the global galactic wind (in Gyrs) - columns (4)-(6): the total 
mass of gas, oxygen, and iron expelled by the system until the present epoch (to = 12 
Gyr) - columns (7)-(9): the present-day photometric properties of interest for this paper 

- column (10): the V-band luminosity- weighted metallicity - column (11): present-day 
mass in luminous matter (star+remnants+gas) - column (12): the total present-day mass 
(luminous plus dark, assuming initial dark to luminous mass ratio of ten) to luminosity ratio 

- column (13): the mass fraction of gas expelled by the system over its lifetime (relative to 
the total initial mass) - column (14): the final system's gas mass fraction (relative to the 
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final luminous mass). 

EDITOR: PLACE TABLE | HERE. 

Recall that the flatter x = 0.95 IMF was chosen for consistency with our earlier work 
Matteucci & Gibson (1995), which used predominantly the code of Matteucci & Tornambe 
(1987), with the Matteucci (1992) extensions. Our "minimal" models in the table above 
are similar to our earlier work, although we note that the later wind times found in our 
earlier work is due mostly to our using the classic Model Ao evolution for the SNe energy, 
as opposed to Model B3 ( |Gibson 1994 bp . Figures [I] and |2| illustrate that our template is 



successful at reproducing the mean of the colour-luminosity-metallicity relations observed 
in present-day ellipticals and dwarfs, over ~16 magnitudes in V. As in Yoshii & Arimoto 
(1987), the astration parameter is treated as a free parameter in order to recover the 
observations. For masses M g J> 10 7 M Q , v oc M~ ai , which is the expected behaviour if the 
initial time-scale for star formation is set by the mean collision time of fragments in the 
collapsing proto-galaxy ( [Arimoto fc Yoshii 19871 ; [Matteucci 1994 ). 



We mention in passing that the two dwarf ellipticals in Figure [l] for which the model 
apparently overestimates the metallicity (NGCs 205 and 147) may themselves only be 
lower limits. Yoshii & Arimoto (1987) speculate that this is due to that fact that they 
were based upon individual red giants in the galaxies outer halos, and as such, may not be 
representative of the overall metallicity of the system, and in particular, the cores of the 
galaxies. 

Our predicted M t /L v , for the standard model, is proportional to Ly 0m , which is 
flatter than the M t /Ly oc £-°- 13 ^-°- 31 observational constraint from the studies of dark 
matter scaling laws ( |Kormendy 199C ). It's important to stress that we have not attempted 



any "fine-tuning" of the input dark matter distribution in order to recover the apparently 
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steeper observational relation. A global input ratio of 10:1, both in mass and radius (dark 
to luminous), naturally led to the -0.07 slope. It is certainly easy to recover the steeper 
value by systematically increasing/ decreasing the dark-to-luminous initial mass ratio as one 
moves to smaller /higher initial masses, but for this work we have avoided any tinkering 
in that direction, especially because we wanted to be in the best situation for having the 
maximum mass ejection from dwarf galaxies and also because of the uncertainties relative 
to the amount of dark matter in galaxies. It is gratifying to note, at least, that the ratio 
does increase with decreasing luminosity, contrary, for example, to what was found in the 
earlier models of Yoshii & Arimoto (1987), which admittedly did not include any dark 
matter halos (see their Figure 6). Finally, our Mi nm /Ly ratios for giant ellipticals are all 
consistent with the observed values of ~ 10 ± 5 (H = 85 km/s/Mpc)- observational values 
of the M/L ratios refer to the internal parts of ellipticals where the dark matter is not 
evident- and increase slightly with Ly [Mi um jLy oc Ly 1 ). This behaviour is observed in 
real elliptical galaxies although the exponent is ~ 0.2 QBender et al. 1992| ). A variation 



of the IMF from galaxy to galaxy could steepen the predicted relation ( [Matteucci 1994| ) 



without any remarkable consequence on the amount of ejected matter from galaxies. On 
the other hand, the trend of the mass to light ratio could be due to an increasing amount 
and/or concentration of dark matter with galactic mass ( [Renzini fc Ciotti 1993|) . 



We re-iterate the point made in Section [2.1| regarding the scatter in the colours for 



the dwarfs shown in Figure [2| Our models trace the lower envelope of the dwarfs. We 
recognise that many of these systems show signs of intermediate age populations ([Ferguson 



fc Binggeli 1994j) , indicative more complex star formation histories than we are capable of 



modelling with the existing version of our code. The lower envelope represents the predicted 
colours for those systems which have the simple single early burst of star formation. Those 
which have subsequent phases of star formation will scatter redward in the plot as despite 
being younger, they will occur in gas that has been pre-enriched by previous generations of 
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dying stars. 

We have not tabulated the [Mg/Fe] in the underlying stellar populations, but we state 
here that the giant ellipticals have a luminosity-weighted [Mg/Fe]~ +0.35, consistent with 
the mean observed values from Worthey et al. (1992). Our [Mg/Fe] values are almost 
constant across the luminosity range covered in their sample, whereas their data implies a 
gentle increase in the value as a function of increasing luminosity, albeit with a great deal 
of scatter. This increase in the magnesium overabundance with respect to iron has been 
addressed recently by Matteucci (1994). 

We have also not shown the fraction of the thermal energy at tow which is due to 
thermalised kinetic energy from pre-SN stellar winds in massive stars versus that due to 
thermalised SNe ejecta. For giant ellipticals this "wind" contribution is ^ 5%, and for low 
mass dwarf spheroidals, ~ 10 — > 20%. This only corresponds to an < 3% reduction in tow 
for giant ellipticals (^ 10% for dwarfs) when compared with models run with no energy 
input from stellar winds. Adjusting v by a few percent from the values listed in column 
(2) of Table [I] would compensate easily for any minor difference in wind epochs (and hence 
the resultant photo-chemical properties) for models run with and without stellar winds, 
further illustrating the point made in Gibson (1994a, 1996c) that stellar winds do not play 
an important role in driving the galactic wind in models of this ilk. 

One last input ingredient we wish to touch upon further is the initial mass function. 
Arimoto & Yoshii (1987), David et al. (1991), and Matteucci & Gibson (1995), amongst 
others, have all put forth persuasive arguments for a "flatter than Salpeter (1955)" IMF 
in elliptical galaxies. Some arguments are based upon ICM abundances (e.g. the latter 
two references), some upon the implications for the underlying ellipticals' photometric 
properties (e.g. the first reference). The situation is discussed further in Gibson (1996a), 
but we wish to at least draw attention to some interesting aspects of the IMF selection. 
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In Table [| we show a comparison of our template x = 0.95 IMF with that of the 
steeper, canonical x = 1.35 IMF of Salpeter (1955), each with the same upper and lower 
mass limits, as before. We only show one initial mass (M g (0) = 10 12 M Q ), and one post-tcw 
model - the "minimal" model - for succinctness. The astration parameter v has once again 
been tuned to reproduce present-day ellipticals which follow the observations of Figures [l] 
and 0. Similarly, the binary parameter A has been chosen a "posteriori to ensure consistency 
with the average present-day SNe la rates mentioned in Section |2.2| (specifically, Ri a = 0.12 
SNu), although the slight difference in the binary parameter A (0.030 versus 0.045) is 
not important to the results. Column (1): the IMF slope, by mass, x - column (2): the 
astration parameter v - column (3): binary mass fraction (3 — > 16 M ) - column (4): the 
galactic wind time - columns (5)- (7): the masses of gas, iron, and its [O/Fe] abundance, 
ejected at tew _ column (8): the mass fraction of gas expelled at tow ~ column (9): [Mg/Fe] 
in the underlying stellar population. 

EDITOR: PLACE TABLE | HERE. 

The key thing to note from Table ^| is that selecting a steeper IMF, but retaining the 
same astration parameter as for a flatter slope, leads to colours which are significantly 
bluer (~ 0.4 magnitudes in V-K for this example), and metallicities which are significantly 
lower (~ 0.3 dex), than those observed. A more "gentle" star formation scenario (the 
third line in the table) is required in order to allow the enrichment to proceed to such a 
level as to match the observations. This results in a factor of ~ 4 decrease in the absolute 
quantity of gas mass and oxygen mass ejected at tew, and a factor of ~ 2 decrease in the 
iron mass ejected. Perhaps, more importantly, the stellar [Mg/Fe] with this "template" 
x = 1.35 IMF is inevitably driven down to the solar ratio, no longer in agreement with the 
Worthey et al. (1992) observations of ~ +0.2 — > +0.3. In a related vein, the [O/Fe] ratio 
in the ejected gas is also only the solar ratio. Since this model is only marginally greater 
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than an galaxy, and we have not included any post-tcw la-driven contribution (which 
could only drive this ratio further downward), it would not bode well for any attempt at 
reproducing the oxygen-to-iron overabundance seen in the x-ray observations of the ICM 
gas. It is supporting evidence such as these last points which leads us to conclude that an 
IMF somewhat biased towards high mass stars (i.e. a slope of x ~ 1) is a necessity in the 
wind models adopted to date, in agreement with our earlier study, which was based upon 
the chemical properties alone ( jMatteucci fc Gibson 199"5| ). 



EDITOR: PLACE FIGURE | HERE. 



EDITOR: PLACE FIGURE § HERE. 
EDITOR: PLACE FIGURE HERE. 

EDITOR: PLACE FIGURE § HERE. 

Using the photo-chemical properties of the ejecta, and the resultant galaxies, shown 
in Table [I], we can use equation ^ to determine the predicted ICM mass of an element 
(or simply the gas mass) originating in the elliptical/spheroidal population of a cluster 
of luminosity Ly^ S0 . This is a particularly nice aspect of the formalism - specifically, the 
predicted ejected mass is normalised to the actual observed cluster luminosity, and does 
not need to be inferred from integrating over the mass functions, assuming some galactic 
mass-to-luminosity relationship ( [Arnaud et al. 1992]) . 

Figure ^ shows the results of said analysis for iron for the minimal (i.e. no post-icw 
contribution to the ejecta) models of Table [I]. The solid curve shows the predicted iron 
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ICM mass as a function of cluster luminosity, assuming the elliptical/lenticular population 
is well-described by a single Schechter (1976) luminosity function with M v = —20.6 and 
a = —1.45. The model is only marginally consistent with the most-metal poor observed 
distribution (shaded region), at a given luminosity. In general, the predicted iron mass is 
~ 1/2 that of the mean of the observed relation. 

Now, let us follow the suggestion of Trentham (1994) (and the supporting observations 
discussed in Section |2.1|) and take the cluster dwarf spheroidal population to be represented 
by a separate Schechter (1976) luminosity function from that of the giant elliptical, with 
My = —16.7 and a = —1.90. We now assume that 30% of the cluster E+SO luminosity 
originates in this dwarf population, which is the upper limit set by diffuse background light 
measurement in rich clusters of galaxies (e.g. [Melnick et al. 1977fe |Thuan fc Kormendy 



1977 ). The remaining 70% is associated with the giant luminosity function, with My and 
a as in the previous paragraph. As the global cluster luminosity does not change, we are 
effectively shifting the available light from one part of the luminosity function to another. 

The top, heavy dotted, curve in Figure | shows the result of this 70/30 distribution. 
The predicted ICM iron mass is now only ~ 76% that of the single luminosity function 
result. The reason for this is as just implied - we have shifted the emphasis, somewhat, 
from giants, which for the minimal model, eject ~ 0.008 M of iron per unit luminosity, to 
dwarfs, which at the low mass end, only eject ~ 1/40 this amount. In this 70/30 luminosity 
split, the dwarfs only contribute ~ 10% the absolute mass of iron that the giants do. 

Figure |] shows the corresponding prediction for gas, in the minimal model. Here we 
see that the single luminosity function can only account for ~ 18% of the observed ICM gas 
at a given cluster luminosity. The double LF is able to raise this to ~ 23%, for much the 
same reason that it led to lower iron - specifically, the ejected gas mass per unit luminosity 
is ~ 5.5 times higher at the low mass end of the minimal models in Table [I], as compared 
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to the giant end. By shifting the luminosity distribution to favour the low mass end of the 
spectrum, we do end up boosting the predicted cluster ICM gas mass which originates, but 
not nearly to the level necessary to explain 100% of the gas. 

Because there is no post-tow contribution to the winds, and the wind epochs are quite 
early (^0.1 Gyr), it is not surprising to note that the predicted ICM gas [O/Fe] is ~ +0.4, 
consistent with a prominent Type II origin to the elements ([Woosley fc Weaver 19951 ), 
and easily consistent with the observed oxygen overabundance relative to iron observed by 
ASCA ( Mushotzky 1994|) . Besides the slight underproduction of iron just alluded to, a 



more notable problem with the minimal model is the predicted final fraction of gas in the 
system at the present-day. From column (14) of Table II], we can see that the ratio of gas 
mass to total luminous mass, at tQ = 12 Gyr, ranges from ~ 0.6 (dwarfs) to ~ 0.4 (giants). 
This is considerably higher than the observed maximum of ~ 10% ([Forman et al. 1985|) . 
Obviously then, some other mechanism has to come into play. We shall return to this point 
momentarily. 

The parallel analysis for the maximal and standard models are shown in Figures [5] 
though [8]. The maximal model overproduces iron by a small amount (~ 50%), but is still 
within the observed maximum for a given cluster luminosity. Because of the increased 
importance of the Type la SNe in the post-tcw ejecta, the [O/Fe] of the ejected gas is 
decreased by ~ 0.3 dex, when compared with the minimal model of Figure |[ The difference 
between the single and double luminosity function scenarios is negligible, although the 
dwarfs do contribute ~ 1/3 of the total iron in the maximal model, as opposed to the 
~ 1/10 we saw in the minimal case. 

Of more interest for the maximal model, perhaps, is the result shown for the total 
gas mass, in Figure ^. Recall that for the minimal model of Figure ^, the ellipticals in 
a cluster could contribute ~ 18 — > 23% of the observed ICM gas. Now, considering the 
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maximal model, we see that for the single/double luminosity function scenario that the 
cluster ellipticals can contribute approximately 33/38%, of which ~ 40% derives from the 
"dwarfs". This is one of the key conclusions of our work - even assuming an extreme 
scenario in which all the gas returned by dying stars is ejected continuously to the ICM, 
neither the giant ellipticals nor the dwarf spheroidals (nor their sum, for that matter) can 
be responsible for all the gas observed in the ICM of galaxy clusters, provided that we insist 
that the resultant galaxies reflect the photo- chemical properties of present-day systems. 

Figures |7| and § show the results for our standard model. The iron mass predicted for 
our model ICMs lies along the median of the observed region, but once again, the gas mass 
accounted for is at most ~ 35%, with a similar distribution of giant/dwarf origin as was 
found for the maximal model. The predicted [O/Fe] of ~ +0.2 is consistent with the ASCA 
observations ( [Mushotzky 1994[) . The dwarfs are identical to the maximal models. This is 
not surprising as their shallow potential wells facilitate the continual expulsion of gas during 
the post-tow regime, a result reflected in the recent hydrodynamical simulations of Wang 
(1995). The giants (i.e. M g (0) = 5 x 10 10 and 1 x 10 12 M ) in this "standard" model only 
drove steady winds for 0.22 and 0.43 Gyrs, respectively, beyond which the binding energy 
of the returning gas from dying stars "catches up" with the more rapidly cooling term for 
the SNe- heated gas. At this point, one might be tempted to re-ignite star formation, and 
indeed this is exactly what Arimoto (1989) and Ferrini & Poggianti (1993) consider. As this 
introduces further free parameters into the picture, we decided not to pursue this option, 
and we simply direct the reader to their excellent studies of late-time star formation within 
the framework of the wind model. 

In this vein though, we note that the final gas fractions for our massive model ellipticals 
range from ~ 20 — > 25%, which is still higher than the < 10% expected from observations 
( [Forman et al. 1985 ). We are not overly concerned by this apparent discrepancy, partly 
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because of the uncertainties in predicting the exact mass of post-tcw ejecta, and partly 
because of the neglection of post-tcw star formation. Indeed, the high final gas fractions 
were similarly a problem with the original models of Arimoto & Yoshii (1987), and as 
shown in Arimoto (1989), recurrent periods of star formation (at a much reduced level 
to that in the initial burst) are a natural consequence of the late-time cooling of the gas, 
although we re-iterate that the exact amount of both the post-tcw star formation and gas 
ejection is highly sensitive to the assumed fraction of thermal energy in the gas immediately 
after tew- Arimoto (1989) alleviated his ~ 20% gas fraction "problem" without unduly 
altering the present-day photo-chemical properties of the galaxies, with limited post-tcw 
star formation/gas ejection. Again, we do not consider such evolution in our models, 
but anticipate that such a scenario, identical to his, would similarly reduce our final gas 
fractions without altering our photo-chemical properties. Still, this is an obvious avenue for 
future research. 

Column (13) of Table [3] shows the initial mass fraction which is ejected into the ICM 
in the form of gas. Regardless of model-type, the values range from 7 « 0.04 — > 0.08. 
These compare favourably with the lowest value considered by Trentham (1994), but are 
not consistent with his higher value of 0.33. It is important though to recall that our values 
were derived from a self-consistent wind model which leads to ellipticals (giant and dwarf) 
with photometric and chemical properties which reflect those observed at the current epoch, 
whereas the values considered by Trentham (1994)were somewhat arbitrary. 



4. Summary 

In summary, we stress that by considering coupled photometric and chemical evolution 
of simple spheroidal models within the framework of a galactic wind model, the precursors 
to the dwarf spheroidals which we observe today in clusters do not seem capable of 
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providing the gas necessary to explain the mass observed in the ICM. Our standard model 
with a giant elliptical luminosity function contributing 70% of the cluster luminosity, and 
the remaining 30% originating in a low luminosity diffuse contribution, shows that the 
precursor dwarfs contribute at most ~ 15% of the ICM gas, with the giants responsible 
for ~ 20%. The remaining ~ 65% is then, by definition, "primordial" in origin. This is 
somewhat less than the ~ 80 — > 90% first proposed by Matteucci & Vettolani (1988), so 
we conclude that the precursors to the dwarfs do contribute a non-negligible amount of 
gas to the ICM of clusters, but that this amount is nowhere near enough to explain all the 
gas. Our self-consistent treatment led to ejected initial mass fractions of 7 ^ 0.08, allowing 
us to exclude the 7 = 0.33 assumption of Trentham (1994). Whether this non-galactic 
wind-originating gas has resided in the cluster since early epochs, or supplied by later gas 
infall, or simply failed to completely collapse into galactic potentials is still a matter under 
vigourous investigation (e.g. White et al. 1993, Nath fc Chiba 1995] , and references therein). 



We admit that the parameter space in these wind models is large, but we explored a 
much larger parameter space than presented in this work and found that without resorting 
to extraordinarily ad hoc assumptions regarding one or more of the input ingredients, we 
can not envision a scenario which will eject the required gas mass and still honour the 
underlying photo-chemical properties of the resultant galaxies and the ICM abundances. 
Moreover, we stress the point that our estimate for the total gas ejected from galaxies is an 
almost maximal one, so that any reasonable change of parameters goes in the direction of 
further decreasing the amount of matter which can be restored from galaxies into the ICM, 
thus reinforcing our conclusion. 
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a Template Models: Minimal models suppress all post-tow SNe la-driven winds; maximal models eject all returning gas 
in said phase to the ICM; standard models include low mass systems which follow the maximal scenario, and high mass 
systems which only drive winds for ~ 0.2 — > 0.5 Gyrs subsequent to tow- Nucleosynthetic yield sources are Woosley & 
Weaver (1995) for high mass stars, Rcnzini & Voli (1981) for single low and intermediate mass stars, and Thielemann, 
Nomoto & Hashimoto (1993) for binary star supernovae Type la. Luminosities and colours derived from Worthey's (1994) 
isochrones. Supernovae remnant thermal energy follows model B3 of Gibson (1994b), which is derived from Cioffi, McKcc 
& Bertschinger (1988) and Larson (1974). A single power law initial mass function, by mass, of slope x = 0.95, and lower 
and upper mass limits of 0.2 M Q and 65.0 M , respectively, was used, v is the astration parameter for the star formation 
rate in equation (1). the galactic wind time tew is m units of Gyrs. A binary parameter A = 0.03 was adopted. Initial 
dark-to-luminous mass and radial extents of 10 were chosen. 7 is the initial mass fraction ejected in the form of gas, and 
/g as is the gaseous-to-luminous mass fraction at tc — 12 Gyr. 
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Table 2. 10 12 M Model - IMF Comparison 1 
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a A comparison of the predicted ejecta composition for two different IMF slopes for an initial gas mass 
of 10 12 M©. The x = 0.95 entry is identical to the template model of Table 1. A Salpeter (1955) x = 1.35 
slope is shown for comparison - one with the identical v as the template, one with v set to recover the 
colours and metallicities shown in Figures 1 and 2. Resultant present-day SNe la rates are 0.11 SNu. See 
text and Table 1 caption for model details. 
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Fig. 1. — Observed metallicity-luminosity relation for dwarf (open circles: Smith (1985)) and 
"normal" (open squares: Sil'chenko 1994 - filled circles: Terlevich et al. 1981) ellipticals. The 
[Mg/Fe]=+0.25, Mg2-Z calibration of Barbuy (1994) has been adopted for the normal/giant 
ellipticals. The solid curve represents our adopted template of models from Section 
H = 85 km/s/Mpc is assumed. 

Fig. 2. — V-K colour-luminosity relation for Virgo (open circles) and Coma (filled circles) 
cluster ellipticals and lenticulars (Bower et al. (1992). A Virgo distance modulus of 
(V — My) = 31.54 is assumed, and a shift of A(V — My) = 3.58 has been applied to 
the Coma sample. Ho = 85 km/s/Mpc is assumed. Data for local dwarfs has been taken 
from Thuan (1985). The solid curve represents our adopted template of models from Section 
22. 

Fig. 3. — Shaded region shows the observed correlation between the non-spiral-originating 
V-band cluster luminosity, and the observed ICM iron mass, after Arnaud (1994). Solid 
curve is the single luminosity function model of slope a = —1.45. Dotted lines are the 
components of the two component luminosity function model - the lower curve is the low 
luminosity dwarf spheroidal component with a = —1.90. The middle one is the normal giant 
spheroidal population with a = —1.45. The heavy dotted curve is their sum. 

Fig. 4. — Shaded region shows the observed correlation between the non-spiral-originating 
V-band cluster luminosity, and the observed gas mass, after Arnaud (1994). Solid curve is 
the single luminosity function model of slope a = —1.45. Dotted lines are the components 
of the two component luminosity function model - the lower curve is the low luminosity 
dwarf spheroidal component with a = —1.90. The middle one is the normal giant spheroidal 
population with a = —1.45. The heavy dotted curve is their sum. 

Fig. 5. — Shaded region shows the observed correlation between the non-spiral-originating 
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V-band cluster luminosity, and the observed iron mass, after Arnaud (1994). Solid curve is 
the single luminosity function model of slope a = —1.45. Dotted lines are the components 
of the two component luminosity function model - the lower curve is the low luminosity 
dwarf spheroidal component with a = —1.90. The middle one is the normal giant spheroidal 
population with a = —1.45. The heavy dotted curve is their sum. 

Fig. 6. — Shaded region shows the observed correlation between the non-spiral-originating 
V-band cluster luminosity, and the observed gas mass, after Arnaud (1994). Solid curve is 
the single luminosity function model of slope a = —1.45. Dotted lines are the components 
of the two component luminosity function model - the lower curve is the low luminosity 
dwarf spheroidal component with a = —1.90. The middle one is the normal giant spheroidal 
population with a = —1.45. The heavy dotted curve is their sum. 

Fig. 7. — Shaded region shows the observed correlation between the non-spiral-originating 
V-band cluster luminosity, and the observed iron mass, after Arnaud (1994). Solid curve is 
the single luminosity function model of slope a = —1.45. Dotted lines are the components 
of the two component luminosity function model - the lower curve is the low luminosity 
dwarf spheroidal component with a = —1.90. The middle one is the normal giant spheroidal 
population with a = —1.45. The heavy dotted curve is their sum. 

Fig. 8. — Shaded region shows the observed correlation between the non-spiral-originating 
V-band cluster luminosity, and the observed gas mass, after Arnaud (1994). Solid curve is 
the single luminosity function model of slope a = —1.45. Dotted lines are the components 
of the two component luminosity function model - the lower curve is the low luminosity 
dwarf spheroidal component with a = —1.90. The middle one is the normal giant spheroidal 
population with a = —1.45. The heavy dotted curve is their sum. 
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